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We discuss transport through a lateral quantum dot in the vicinity of a singlet-triplet spin transi- 
tion in its ground state. Extracting the scattering phase shifts from the numerical renormalization 
group spectra, we determine the linear conductance at zero temperature as a function of a Zeeman 
field and the splitting of the singlet and triplet states. We find reduced low-energy transport, and a 
non-monotonic magnetic field dependence both in the singlet and the triplet regime. For a generic 
set of dot parameters and no Zeeman splitting, the singlet-triplet transition may be identified with 
the conductance maximum. The conductance is least sensitive to the magnetic field in the region 
of the transition, where it decreases upon application of a magnetic field. Our results are in good 
agreement with recent experimental data. 



I. INTRODUCTION 

Recent observations of the Kondo effect in semicon- 
ductor quantum dotsi* induced enormous theoretical 
and experimental activity. Compared to magnetic im- 
purities in solids*, quantum dots* have the great ad- 
vantage of tunable couplings, and as a result, in these 
systems the full parameter regime of the single-level An- 
derson impurity model can be systematically explored. 
In addition, new types of Kondo systems have been re- 
alized, like Aharonov-Bohm rings containing quantum 
dots* and multi-level dots*. In vertical dots, a qualita- 
tively new type of Kondo effect associated with a singlet- 
triplet degeneracy has been observed*, which was later 
also explained theoretically 8,9,10 . 

For lateral dots, this singlet-triplet Kondo effect has 
been found as welU*, but with qualitatively new behav- 
ior at low energies, where both the linear conductance 
and dl/dV were found to be non-monotonic. Theoret- 
ical understanding of this behavior has been obtained 
in two limiting cases: Hofstetter and Schoellei— have 
analyzed the interplay between singlet and triplet con- 
figurations with symmetric coupling to the leads, while 
Pustilnik and Glazmani* have considered general asym- 
metric couplings but took only the triplet configuration 
on the dot into account. 

Agreement with the singlet-triplet scenario of 01 nas 
recently been found in another experiment*^, where - in- 
stead of the magnetic field as in - the Stark effect due 
to the gate potential was used to tune the level splitting. 
Another way to control the level spacing and to study 
the triplet-singlet transition by a symmetric (triangular) 
arrangement has been proposed in [l5| . 

Our goal in this work is to extend the studies of [Tiil| 
over the entire range of parameters in a non-perturbative 
way, and to include singlet /triplet degeneracy as well as 
asymmetry of the tunnel couplings in our transport calcu- 
lation. To this purpose we shall use the numerical renor- 
malization group (NRG) in a rather unconventional way: 
We shall extract the T = scattering phase shifts directly 
from the NRG spectrum**** and use these to compute 



the conductance by applying the Landauer-Buttiker for- 
mula. Our calculations are in good agreement with ex- 
periments and reproduce the results in special limits as 
obtained in 02 and [13. 

II. THE MODEL 

In this paper we focus our attention to the lateral quan- 
tum dot system shown in Fig. ^ Such a quantum dot is 
usually formed by gate depletion in a two-dimensional 
electron structure. The conductances between the dot 
and the leads can be controlled by applying a voltage on 
the gates separating the dot region from the leads. In this 
system Coulomb blockade develops when the dimension- 
less conductances gh,R = Gl.r/ '(2e 2 '/h) between the dot 
and the left and right leads drop below 1: In this regime, 
electrons have to pay typically an energy of the order 
of the charging energy ~ Ec to get onto the dot, and 
therefore transport through the structure is suppressed 
when the temperature drops below Ec- Since Coulomb 
blockade develops when the lowest propagating mode in 
the point contacts (between dot and leads) is at pinch- 
off, all higher modes can be neglected. The leads can 
thus be modeled by a single conduction electron mode 
each**. Note that this is specific to lateral dot systems: 
In vertical dots (see e.g. 001 X many conducting modes 
participate in transport. 

The Kondo effect may appear at temperatures T <C 
Ec if the ground state of the isolated dot is degenerate. 
In this case quantum fluctuations to the leads are rele- 
vant, and may lead to the formation of a Kondo state at 
an energy scale Tk- This condition is obviously fulfilled 
if there is an odd number of electrons on the dot, since 
then the ground state of the dot has a Kramers degen- 
eracy. The Kondo effect then emerges if we lower the 
temperature below the Kondo temperature*?- 

T K ~ SE ex V {-E c /{g L + 9r)5E} (1) 

where SE is the mean single-particle level spacing in the 
dot. 
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Singlet state ( S=0 ) 



Triplet state ( S=l ) 



FIG. 1: Lateral quantum dot geometry. The dot is formed by 
depletion of the two-dimensional electron gas underneath the 
gates. Dark areas denote the positions of the gate electrodes, 
while the black lines indicate the border between depleted 
(white) and undepleted (grey) regions of the 2DEG. 



However, Kondo effect can also be observed in dots 
with an even number of electrons. To understand the 
physics that leads to the Kondo effect in this case, let us 
first study an isolated quantum dot. A generic model of 
an isolated quantum dot can be written as2S 



H, 



dot 
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s is the total number of electrons 
£w d ns (*ss>/2) d„., is the to- 



where N = £„ s d{ s d r , 
on the dot, and S = 
tal electronic spin on the dot. The operators d^g create 
electrons on a single particle level of the dot, labelled by 
the spin s and a discrete quantum number n. Eq. J2J de- 
scribes the electron-electron interaction at the mean field 
level. In general, more complicated interaction terms 
should be present in H ao ^. These terms are, however, 
relatively small for dots with a large number of electrons 
and furthermore, they do not influence our discussion of 
the singlet triplet transition below, therefore we shall ne- 
glect them. 

The parameter N in Eq. denotes the dimensionless 
gate voltage, and it sets the average number (N) of elec- 
trons on the dot, B is the Zeeman field, and Ec and Es 
stand for the charging energy and the Hund's rule cou- 
pling, respectively. In the rest of the paper we shall focus 
our attention to the case where Nq is close to an even inte- 
ger, and we shall assume Ec 3> 8E > Es, characteristic 
of lateral dots with a large number of electrons^. Under 
these conditions, the ground state of Eq. (|2J) is typically 
a singlet state shown on the left in Fig. |3 and therefore 
no Kondo effect occurs. If, however, the last occupied 
(e_i) and the first empty state (e+i) happen to be close 
enough to each other, so that A = e+i — e_i < 2Es, then 
the system will form a triplet state to gain energy from 
the Hund's rule coupling by rearranging the level occu- 
pancy (see Fig. |2J) . In this case the ground state is three- 
fold degenerate and a Kondo state can be formed. This 
transition occurs when the energy difference between the 
singlet and triplet states approximately vanishes: 

A S t = e+i - £-i - 2E S « . 
The Kondo effect is driven by virtual charge fluctua- 
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FIG. 2: Possible triplet and singlet ground states of an iso- 
lated dot described by Eq. 



tions to the leads. Since the coupling between the leads 
and the dot is weak, gx, gn -C 1, we shall describe these 
fluctuations in terms of the following simple tunneling 
Hamiltonian: 

Pleads = J2tk a lks a a ks' u = R,L (3) 

aks 

■^tunneling = ^an a aks^ns + H.C. (4) 



an ks 



The operators a aks in Eq. @ create electrons with mo- 
mentum k, energy and spin s in the left and right 
leads (a = L,R), respectively. Eq. (0} describes the tun- 
neling between the dot states (labelled by the spin s and 
a discrete quantum number n) and the leads. 

In the vicinity of the singlet-triplet transition charge 
(and spin) fluctuations on the dot are dominated by fluc- 
tuations to the states n = ±1. It is therefore sufficient 
to restrict the summations in Eqs. (0) and J2J and keep 
only these two states: 



n A/2, 



±1. 



(5) 



This projection will not crucially influence our results, 
and the two-level model is expected to capture all uni- 
versal aspects of the problem until the relevant energy 
scales of the system that govern the singlet-triplet tran- 
sition (i.e., the Kondo scales, B, and Ast) are all small 
compared to all other energy scales (Ec, SE etc.) in the 
problem. However, when projecting to these two states, 
one also has to keep in mind that at energy scales larger 
than max{<5-E, Ec} fluctuations the the excited states of 
the dot destroy those coherent contributions that give 
rise the Kondo effect. Therefore, at the same time the 
effective band-width in Eq. Q has to be reduced to 
D = max{8E,E c )22.. 

In the following sections we shall study the conduc- 
tance of the dot when one gradually drives the system 
through the triplet singlet transition. There are several 
experimental techniques to achieve this goal and tune A 



3 



and thus Ast- The easiest method is to apply a magnetic 
field H± perpendicular to the plane of the dot. The pri- 
mary effect of H± is to change the orbital wave functions 
and thus to control the value of A(iJj_)£Ai By lineariz- 
ing A(H±) in the vicinity of the transition, one can thus 
attempt to make a direct comparison of the experimen- 
tal data with the (calculated) linear conductance across 
the system G — (2e 2 /h)g(A). Unfortunately, H± also 
induces a Zeeman field B(H±): While this Zeeman split- 
ting B remains small compared to the orbital effect be- 
cause of the smallness of the electron effective mass and 
the <?-factor in GaAs 2 ^, yet it is comparable to the exper- 
imentally observed Kondo scale on the triplet sideii, and 
leads to some complications when trying to make a direct 
comparison with the experimental data. The difficulties 
related to Zeeman splitting were avoided in the experi- 
ment 0| by tuning the singlet-triplet transition via the 
Stark effect due to an inhomogeneous gate potential. 

Another way to avoid this problem is to use triangular 
dots where the degeneracy of the states e±\ can be di- 
rectly controlled by changing the shape of the doti^. To 
investigate the dependence of <?(A, B) on B experimen- 
tally, one may in addition apply a strong in-plane field 
Hit , which only generates a Zeeman splitting without in- 
fluencing A essentially. 

With the method outlined in Section IIVI we cannot 
compute the finite temperature conductance of the dot. 
However, the temperature (T), source-drain bias (V), 
and Zeeman field (B) dependences of the conductance are 
expected to be qualitatively (but not necessarily quanti- 
tatively) similar to each other. Therefore many of our 
results can be used to understand qualitatively the finite 
temperature and finite bias behavior of the dot in the 
absence of the Zeeman field B by replacing B with V or 
T. 



III. LINEAR CONDUCTANCE 

As discussed above, in the vicinity of the singlet-triplet 
transition the Hamiltonian of the dot can be truncated 
to that of a two-level system (Eqs. (0) and (0). The 
tunneling, Eq. Q), couples the two levels to the two leads, 
and the four tunneling amplitudes t a%n form a 2 x 2- 
matrix 

t= ( t L > +i ;*» +i Y 

Consider first the special case when one of the eigenval- 
ues of the matrix t is zero, while another one is finite^. 
Obviously, in this case the dot effectively interacts only 
with a single species (a single "channel" ) of conduction 
electrons. When the dot is in the triplet state, a sin- 
gle electronic channel can screen only half of its spin 24 . 
Accordingly, the system should exhibit a quantum phase 
transition 25 : the ground state changes its symmetry from 
a singlet to a doublet as A decreases below a certain 
critical value Ac ~ E$- The linear conductance g(A,0) 



across the system at T = and B = is then strongly 
A-dependenO 2 .: g(A,0) cx 6(A C - A). At A < A c 
(when the dot is in the triplet state) the conductance is 
a monotonically decreasing function of B. To the con- 
trary, at A > Ac the conductance as a function of B 
first increases, and then drops with the increase of B. At 
a fixed finite B the conductance <?(A, B) is expected to 
be a smooth non-monotonic function of A with a broad 
asymmetric peak near the transition. 

8 A T=B=0 
triplet singlet 




/ A 
A = 2E S 

max s 

FIG. 3: Sketch of the zero temperature conductance as a 
function of the splitting between the two dot levels for generic 
tunneling amplitudes. 

In a general situation, however, det i ^ 0, and both 
eigenvalues of the matrix t are finite. In this case, the 
dot is coupled to two electronic channels, which is suffi- 
cient in order to completely screen the dot's spin 2 ^ even 
in its triplet state. As a result, the ground state of the 
system is a singlet at all values of A&iS* 2 ^. Accordingly, 
when the dot is coupled to the leads, the singlet-triplet 
transition turns to a crossover. As we discuss below, 
the conductance g(A, 0) at T = and B = is then 
expected to be a smooth non-monotonic function of A, 
slowly (logarithmically) decreasing-^ at large |A — 2E S \ 
(see Fig. 121. Far away from the crossover, the conduc- 
tance is a non-monotonic function of B: it first raises 
and may reach a value g ~ 1, and then drops again as 
B increases. This has been shown for both the triplet 
side of the crossoveni^Si (A <C Es), and for the singlet 
(A > E s ) side of For a fixed small B, how- 

ever, the conductance g(A, B) can be an even more com- 
plicated function of A with a local minimum at some 
A < Es and maxima at the triplet side of the transition 
(A < E s ) and around A ~ E s . 

In this paper, we focus on a particular limit when 
both eigenvalues of t are finite, but still very different: 
| det 1 1 <C (Tr t) . In addition, we constrain our consid- 
erations to a special subse1i2§*21 of the matrices t that can 
be parametrized as 

i=±( V+1 V+1 ), WO. (6) 

This tunneling matrix arises, e.g., for a fully symmet- 
ric arrangement when the states n — ±1 are symmet- 
ric/antisymmetric under reflection. The choice 10, al- 
though by no means general, nevertheless captures the 
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essential physics of the system. At the same time, as we 
demonstrate shortly below, the choice © allows one to 
express the conductance g(A, B) at T = in terms of 
simple phase shifts that can be extracted from the finite 
size NRG spectra^, or in some cases can be related to 
thermodynamic quantities^Sl. 

Since the ground state of the system is not degenerate, 
electrons scatter elastically at T = 0. The scattering am- 
plitudes S s:a . a > of an electron with spin s from lead a' to 
lead a form the scattering matrix S sa - S ' a ' = S ss > S S ; a ,a'- 
The 4x4 unitary matrix S can be diagonalized by a ro- 
tation in the R — L space to the new basis of channels 
n = +1 and n = — 1, 



USU* = diag{e 2i5 " s } , U = e 



(7) 



Here t % are the Pauli matrices acting in the R — L space. 
The linear conductance at T = is related to the off- 
diagonal elements of S by the Landauer formula, 



c2 I 



and, by making use of Eqs. can now be expressed via 
the scattering phase shifts 5 ns at the Fermi energy. In 
general, the angle i? and thus g depend explicitly on the 
parameters of the microscopic Hamiltonian. However, 
with the choice Eq. I© the tunneling part of the Hamil- 
tonian can be trivially diagonalized with the introduction 
of even and odd states (corresponding to # = 7r/4, tp = 0) 



C+l,ks 
C-l ks 



_ _J_ ( aRks 
y/2 V aRks 



and the Hamiltonian takes the form 



H = 



nks C nks + V n [ C nks^ns 



dLks 
O-Lks 



H.C. 



H. 



dot 



(8) 

The Hamiltonian (jSJ with H dot given by Eqs. ©,10 ob- 
viously conserves the total number of fermions in each 



channel, Af n = ^ A 



^nks ^nks 



- $ d 



implying that 



the scattering matrix is diagonal in n. 

In this case the conductance simply reads 13,29 : 



9 = 5 E sin2 ^+ 1,s ~ 5 - 1 -^ ■ 



(9) 



Following Nozieres, at small B the phase shifts can be 
expanded to linear orde»2fl^i. As explained in |30j , a 
simple analysis shows that in the presence of a Zeeman 
field the conductance at low fields behaves as 



«/(A,S) w 5 (A J Q) + [l-2 5 (A > 0)] 



B 



B*(A) 

where the scale B* is a A-dependent energy scale and 
. 9 (A,0)=sin 2 p +1 (A)-MA)]. 



(10) 



In other words, far away from the transition the con- 
ductance should increase with B for g(A,0) < g* w 1/2 
while it should typically decrease in the transition region, 
where g(A, 0) > g*. We have to emphasize that the value 
g* ps 1/2, see Eq. (fTUl) . is specific to the choice ©, and 
not universal. 

In the following, we shall compute the phase shifts us- 
ing the powerful numerical renormalization group and 
determine the conductance using Eq. @. 



IV. NUMERICAL RG, FINITE SIZE SPECTRA, 
AND PHASE SHIFTS 



The Numerical Renormalization Group (NRG) was 
originally developed by Wilson to solve the Kondo 
problem 32 . In contrast to scaling approaches, it is non- 
perturbative and does not encounter any logarithmic sin- 
gularities at low energy. Since then it has been extended 
to the calculation of dynamical quantities 3 ^ 3 ^ and has 
been applied to a variety of quantum-impurity problems, 
including transport calculations for interacting quantum 
dotsiL^. 

In principle, the NRG method enables us to compute 
the expectation value of any local operator. Thus one 
plausible method to obtain the phase shifts would be to 
compute the occupation numbers N ns — (<ijj s d ns ) in the 
ground state by NRG and then use the Friedel sum rule 
to get the phase shifts as&Si: 



S ns = ttN„ 



(11) 



Unfortunately, this procedure turns out to deliver un- 
physical results, especially on the triplet side. The origin 
of this problem is probably related to the relatively small 
bandwidth in our calculations and the large logarithmic 
corrections that appear on the triplet side: As is obvi- 
ous from the derivation of the Friedel sum rule^, there 
is, in general, a correction term to Eq. (|ll|l which van- 
ishes in the infinite band width limit. In our calculations, 
the Kondo temperature is comparable to the bandwidth 
D. Furthermore, on the triplet side of the transition the 
convergence to the ground state is rather slow due to the 
large ferromagnetic residual exchange coupling, generic 
to underscreened Kondo models, that can even lead to 
singular non- Fermi liquid properties 3 ^. Therefore we can- 
not use the Friedel sum rule to determine the conductiv- 
ity of the dot. Nevertheless, the Friedel sum rule should 
definitely work in the limit of very large cut-offs. It there- 
fore gives us a useful tool to understand qualitatively the 
evolution of the phase shifts and the conductance. 

Fortunately, as discussed below, we can determine the 
phase shifts directly from the NRG spectra without mak- 
ing use of the Friedel sum rulei&. This method has been 
first used to compute the conductance of a mesoscopic 
double dot systemii. To do this, we only have to know 
that the model flows to a universal Fermi liquid fixed 
point that can be characterized by four phase shifts. 
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To explain the method, let us briefly discuss Wilson's 
NRG procedure. The key idea of Wilson is to map the 
original problem to a semi-infinite chain. After this 
transformation our two-channel Hamiltonian maps onto 
a system of two parallel tight -binding chains coupled to 
the dot 

H = H dot + v n (dl s c Q ns + h.c.) 

71,8 

oo 

where the hopping matrix elements decay exponentially, 
^ sa A±A — A -1 / 2 , with A the NRG discretization param- 
eter. In our calculation we take A = 3. The Ci <ns rep- 
resent conduction electron excitations at a length scale 
~ A 4 / 2 . The NRG assumes a flat dispersionless density of 
states for the electrons and all energies are measured in 
units of the bandwidth D. The operator co. ns is simply 
the local field operator defined as co, ns = J k dkckns , 
and it is only this operator that tunnel-couples to the 
interacting impurity part, H dot . Note also that the dis- 
cretization procedure somewhat renormalizes the original 
parameters of the Hamiltonian, which are therefore not 
identical to the ones in Eq. (JBJ. 

Eq. 1121 can be rewritten in a more inspiring way as 

H N+1 = ^ ,2 H N + Y J (^ Nns c N+lns + h.c) (13) 



H = lim 1±A1a-(»-Wh n 

JV->oo 2 



(14) 



where we have introduced the rescaled Hamiltonians 
2 



. A (N-1)/2 H 



N 



Hn = 1 + A-i' 
H N ee H dot + Y v n (di s c ns + h.c. 

71, S 

JV-1 

+ & { C { n s C i+l,na + h - c - 

n,s i—0 



(15) 



(16) 



and Hq denotes the impurity Hamiltonian, H dot , tunnel- 
coupled to co.p S . Eq. (|T4*jl can be solved by iterative di- 
agonalization, by keeping in each step only the lowest, 
most relevant levels. In our calculation, conservation of 
total charge and the z-component of the total spin 



n,s i—0 



n,s i—0 



1/2) 



(17) 
(18) 



have been exploited to increase numerical efficiency. 

A typical example of the spectrum of Hn vs. iteration 
number is shown in Fig. 0] At large iteration numbers 



(i.e. small energy scales), the spectrum remains invari- 
ant, indicating the presence of a stable field theory (fixed 
point) which describes the low energy behavior of our 
model. 

To interpret the spectrum, we have to keep in mind 
that in the NRG procedure the length scale L of the 
system corresponds to L ~ A N / 2 , and therefore, apart 
from some multiplicative factor, the spectrum of Hn is 
the spectrum of a system of size L in units of 2-kvf/L. 
This fixed point NRG spectrum is universal and can be 
characterized by boundary conformal field theories^. In 
our case the fixed point spectrum is that of a Fermi liquid, 
and is therefore identical to that of the following simple 
model 



H 



L/2 



ii, - ^ / dx ipl s (x)(-id x + U ns 5(x))ip ns (x) (19) 



where the ^ ns (x)'s denote chiral fcrmionic fields obeying 
periodic and anti-periodic boundary conditions for even 
and odd A, respectively. This statement means that the 
full many body-spectrum of our model is identical to that 
of the following simple effective Hamiltonian (in Fourier 
space) : 



H 



ip 



Sn,s 



(20) 



where the phase shifts <5„ iS are related to the couplings 
U ns and q takes integer (half-integer) values for even 
(odd) iterations. In more precise terms: 



lim H N = C —Ht p , 

N—>oo 2n 



(21) 



where C is a A-dcpcndcnt constant to be determined 
numerically. The quantum numbers Q and S z in the 




FIG. 4: A typical NRG level flow on the triplet side for 
E c = 1, E s = 0.25, A = 0.2, v +1 = 0.67, «_j = 0.37 and 
B = 0.01. The results were obtained with A = 3 and keeping 
N = 1000 levels. 
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TABLE I: Comparison of the even iteration fixed point NRG 
spectrum in Fig.^Jand that of the effective field theory, Eg. 1201 
We used C — 1.6023, and the phase shifts were — = 0.1757, 

^ = 0.3661, ^ = 0.5843, and ^ = 0.8712. H fp reproduces 
all the quantum numbers correctly and the NRG spectrum 
with a four digit precision. 

fixed point theory are defined as 

Q fP = E : 4,ns%,ns : , (22) 
q,n,s 

Sl P = E S 4,nJi,n S , (23) 
q,n,s 

where : . . . : denotes normal ordering with respect to the 
ground state. 

To illustrate how well this procedure works we have 
enumerated in Table [I] the lowest eigenstates and their 
quantum numbers in the NRG spectrum together with 
those of the effective Hamiltonian Hf p : All quantum 
numbers are in agreement and the eigenvalues match up 
to four digits. In our calculations we have not kept track 
of the conserved quantum numbers A^„. Therefore, while 
we could tell the values of the two spin up and spin down 
phase shifts, S s and S s , in some cases we were unable 
to tell to which channel they belonged. Since, however, 
Eq. 10 is symmetrical under the exchange of n = ±, this 
caused no problem in the calculation of the conductance. 

This method enables us to extract the phase shifts di- 
rectly from the finite size spectrum. The great advantage 
of this procedure is that it works directly on the fixed 
point Hamiltonian and therefore eliminates the problem 
of high energy degrees of freedom that occurred when us- 
ing the Friedel sum rule. This method can be applied to 
any problem where the ground state is that of a Fermi 
liquid. 

V. RESULTS 

In Fig. [3 we have plotted the typical splitting depen- 
dence of the phase shifts for B = and the correspond- 
ing conductance. The qualitative behavior of the phase 



shifts can be understood in terms of the Friedel sum rule 
Eq. Ijlll) as follows: On the triplet side of the transition 
the spin of the dot is completely screened by the electrons 
in the leads, which corresponds to a phase shift tt/2 in 
both channels. On the singlet side, however, both dot 
electrons occupy the level n = — 1. Therefore the phase 
shift in channel n = — 1 should approach (5_i. s ~ n on 
this side, while the phase shift <5+i. s must go to zero. By 
Eq. © this implies that the conductance must approach 
on both sides of the transition while it has a maxi- 
mum <7max ~ 1 around <5_i, s — 5i }S ~ tt/2. Note that the 
transition point is shifted with respect to the bare value 
A = 2Es = 0.5 for the isolated dot due to correlation 
effects. 
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FIG. 5: Phase shifts (top) and corresponding conductance 
(bottom) as a function of orbital splitting A in the absence 
of the Zeeman field at temperature T = 0. Parameters are 
chosen as Ec = 1, Es = 0.25, w+i = 0.67 and v-i = 0.50. 

Application of a Zeeman field B has a more complex 
effect, which to some degree has been discussed in [l3j . 
On the triplet side, A < £ s we have two exponentially 
separated Kondo scales, T£ ^> T£. For Zeeman fields 
B <C T^,T^ all four phase shifts approximately equal 
7r/2, corresponding to a reduced conductance. However, 
for intermediate fields < B < T%, the Kondo ef- 
fect in channel n = — 1 is suppressed, and therefore the 
phase shifts in this channel are expected to be around 
S n =-i,i ~ tv and £ n= _i4 « 0, while 5 n =+i,t ~ 5 n=+ ^ f» 
7r/2, leading to a conductance g ~ 1. Finally, for even 
larger fields, B ^> T^,T^, both Kondo effects are sup- 
pressed, the phase shifts become 8 n =-i,i ~ #n=i,f ~ 
and 5 n —-i^ ~ <5n=i,| ~ tt, and the conductance de- 
creases again to 0. This behavior is shown in Fig. EI 

As we already mentioned in Section \n\ in the exper- 
iment 01 the Zeeman field was non-negligible. It is 
therefore instructive to study the splitting dependence 
of the conductance in the presence of a finite Zeeman 
field B. This is shown in Fig. [7| In a finite magnetic 
field g(A,B — fixed) can show a very surprising A de- 
pendence for some values of B, and may even have a dip 
on the triplet side. As we explain below, this dip ap- 
pears because Tj£ is very sensitive to the distance from 
the transition point, and therefore we may cross over 
from the regime T>, T< > B to T> > B > T< as we 
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FIG. 6: Phase shifts (top) and corresponding conductance 
(bottom) as a function of Zeeman field B on the triplet side 
of the transition at zero temperature. Parameters are Ec = 1, 
E s = 0.25, A = 0.4, v+i = 0.67 and v-i = 0.25. 

are getting farther away from the transition region. For 
intermediate Zeeman fields B the conductance has the 
shape of a smeared step function that becomes sharper 
as the ratio w+ /«_ increases and resembles the one found 
in [12. 




0.2 0.4 0.6 0.8 

A 



FIG. 7: Conductance as a function of A for a fixed value B = 
10 -3 of the Zeeman field and coupling asymmetry v+i/v—i = 
2.2 (top) and v+i/v—i = 8.4 (bottom). The other parameters 
are E c = 1, E s = 0.25, v+i = 0.67. 

The effect of a Zeeman field on the conductance should 
qualitatively resemble that of a finite source-drain bias 
or temperature. We can thus use the -B-dependence of 
the conductance for fixed fixed values of A to identify 
the second energy scale that appears as the width of a 
dip in the zero-bias anomaly in the experiment [TT| . The 
-B-dependent conductance for three different values of A 
is shown in Fig. [S] The conductance first increases for 
small Zeeman fields B on both sides of the transition, 
but it is suppressed for larger values of B. However, it 
decreases monotonically in the transition region, where 
g > g* (see Ea.HU)l. 

We can then identify the low-energy scale B* by fix- 
ing A and measuring the curvature of g(B,A). This is 



FIG. 8: Linear conductance versus in-plane magnetic field 
B for different values of the orbital splitting and the same 
parameters as in Fig. [5] Notice the non-monotonic be- 
havior ("dip") on both the triplet (A = 0.3) and singlet 
(A = 0.6) sides of the cross-over. In the middle of the tran- 
sition (A = 0.52) the conductance is g = 1 at B = 0, and 
decreases monotonically as a function of magnetic field. 



shown in Fig. [5J On the singlet side of the transition 
B* increases approximately linearly with A, and prac- 
tically corresponds to the energy splitting of the singlet 
and triplet multiplcts. Intuitively this can be understood 
in the following way: Applying a Zeeman field one grad- 
ually pulls down the S z = — 1 state of the excited triplet 
state. This becomes degenerate with the singlet state 
if B is approximately equal to the splitting between the 
triplet and the singlet states of the dot. At this special 
value of B a Kondo effect can occur between the ground 
state singlet and the S z = — 1 triplet stated, giving a 
conductance maximum as a function of B. Indeed, we 
can see in Fig. [S] that on the singlet side (A = 0.6) the 
maximum conductance is indeed approximately g ~ 1, 
indicative of the Kondo resonance. Therefore, on the 
singlet side B* must be proportional to the renormalized 
(by quantum fluctuations) singlet-triplet splitting. 

The origin of the dip on the triplet side is very different, 
as explained above. Far away from the transition, the 
spin S of the dot is screened by two consecutive Kondo 
effects. An intermediate magnetic field suppresses the 
Kondo effect with the smaller Kondo temperature, T^, 
and leads to a large, almost unit conductance. There- 
fore, in this region, the energy scale B* should practically 
correspond to the smallest of these two Kondo temper- 
atures, i.e. the one with smaller hybridization v < . Far 
away from the triplet side, TS does not depend strongly 
on the presence of the singlet. However, in the vicinity 
of the singlet-triplet transition, quantum fluctuations to 
the singlet state are extremely important: They mix all 
four dot states, renormalize the value of i><, and gener- 
ate a strongly correlated state with a Kondo scale that 
is determined by the hybridization to the more strongly 
coupled channel, v > . As a consequence, in the vicinity of 
the transition these fluctuations strongly renormalize the 
value of v < and thus influence B* ~ . This is shown 



8 



in Fig. OH 

For the parameters of the figure the second Kondo 
scale would be very small in the absence of fluctuations 
to the singlet state. Therefore, B* is extremely sensi- 
tive to the distance from the triplet-singlet transition, 
which sets the renormalized value of w< , and thus that 
of T£ ~ B*. Far from the transition the effect of the 
fluctuations to the singlet state is small, and they do not 
renormalize TS ~ B* substantially. Thus TS saturates 
on the triplet side and shows only a weak A-dependcncc 
there, in complete agreement with the experimental re- 
sults of Ref. [llj. However, for larger ratios of v+i/v-i, 
the saturation value T^(Ast "C 0) becomes so small that 
it can hardly be observed. 
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FIG. 9: Low-energy scale B* for Ec = 1, Es = 0.25 and 
v+i — 0.67, as determined from the curvature of g(A,B) for 
small Zeeman field B. Note the saturation of B* on the triplet 
side for the larger v~i. 

In Fig. we show the A-dependence of the conduc- 
tance for three different values of B. The conductance 
has a similar behavior for the two lower values of B. In 
these cases, far away on the triplet side, we have a sit- 
uation where T% < B < T^. The smaller Kondo scale, 
Tjf , gradually increases as we approach the transition 
point, and therefore the conductance decreases be- 
comes larger than B, thus giving a dip on the triplet side 
of the transition. 

The third case is very different. Here B — 0.1 is larger 
than both Kondo scales, and therefore the bump on the 
triplet side is suppressed. Even the maximum associated 
with the triplet-singlet transition is shifted to larger val- 
ues of A. We suspect that in this case the bump has an 
entirely different origin and is due to the Kondo effect 
associated with the degeneracy of the S z = — 1 triplet 
state and the singlet stated. 

VI. CONCLUSIONS 

We have studied transport through a lateral quantum 
dot in the vicinity of a singlet-triplet transition in the 
ground state. Our calculation is based on the numeri- 
cal renormalization group, which enables us to calculate 
the full many-body spectrum of the dot. To determine 
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FIG. 10: Linear conductance versus level splitting for differ- 
ent in-plane magnetic fields B and a fixed ratio v+i/v-i = 
2.6. The other parameters are the same as in Fig. |B] 



the conductance we used a new method, first outlined 
in 17], where we extract the conduction electron phase 
shifts with high accuracy directly from the NRG spec- 
trum and use them in combination with the Landauer for- 
mula to determine the conductance at zero temperature. 
The only requirement for this technique to work is that 
the system considered must have a Fermi liquid ground 
state with the phase shifts being well-defined quantities. 

In our calculations have modeled the quantum dot by 
a two-level Anderson impurity Hamiltonian with asym- 
metric couplings to the left and right single-mode lead. 
Our results show that at T = the sharp quantum phase 
transition that occurs for symmetric couplings becomes 
a crossover at finite asymmetry. As a function of level 
splitting, the conductance shows a characteristic maxi- 
mum ( "bump" ) at the crossover point and decreases on 
both the singlet and the triplet side. 

In order to model the effect of finite temperature or 
bias voltage, we have applied an in-plane magnetic field 
B leading to a Zeeman-split triplet. We find that on 
both the singlet and the triplet side, the conductance 
has a non-monotonic behavior as a function of the Zee- 
man field, thus leading to a characteristic "dip" at low 
B. Because of the qualitative equivalence between B and 
finite bias voltage, we expect this structure to appear in 
the differential conductance dl/dV as well. The associ- 
ated low-energy scale (dip width) increases linearly on 
the singlet side, where it corresponds to the renormal- 
ized singlet-triplet splitting. On the triplet side, it de- 
creases rapidly with increasing distance from the degen- 
eracy point, which is in agreement with two-stage Kondo 
screening of the triplet. 

These findings are consistent with the experimental re- 
sults of van der Wiel et al. where a dip in dl/dV was 
indeed found on both sides of the singlet-triplet transi- 
tion. On the other hand, in the measurements of Kogan 
et al. 14], the dip structure was only observed on the 
singlet side, indicating a weak asymmetry below experi- 
mental resolution. 

We expect our results to be relevant also for trans- 
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port through other multi-level Kondo systems like car- 
bon nanotubes or molecules. 
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